Regularizing the r-mode problem for nonbarotropic relativistic stars 
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CNJ . We present results for r-modes of relativistic nonbarotropic stars. We show that the main differ- 

ential equation, which is formally singular at lowest order in the slow-rotation expansion, can be 
regularized if one considers the initial value problem rather than the normal mode problem. flow- 
ed ■ ever, a more physically motivated way to regularize the problem is to include higher order terms. 

This allows us to develop a practical approach for solving the problem and we provide results that 
support earlier conclusions obtained for uniform density stars. In particular, we show that there 
will exist a single r-mode for each permissible combination of I and m. We discuss these results 
and provide some caveats regarding their usefulness for estimates of gravitational-radiation reac- 
' tion timescales. The close connection between the seemingly singular relativistic r-mode problem 

and issues arising because of the presence of corotation points in differentially rotating stars is also 
1 clarified. 
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I. INTRODUCTION 

In the last few years the instability associated with the r-modes of a rotating neutron star has emerged as a plausible 
source for detectable gravitational waves. This possibility has inspired a considerable amount of work on gravitational- 
^Jy wave driven instabilities in rotating stars and our understanding continues to be improved as many of the relevant 
issues are intensely scrutinized (see 0, 13- 0, 0, E| f° r detailed reviews and important caveats on the subject). To 
date, most models for the unstable r-modes are based on Newtonian calculations and the effect of the instability on 
the spin rate of the star is estimated from post-Newtonian theory. This may seem peculiar given that the instability is 
a truly relativistic phenomenon (its driving mechanism is gravitational radiation reaction); but a complete relativistic 
calculation of the oscillation modes of a rapidly rotating stellar model (including the damping/growth rate due to 
gravitational- wave emission) is still outstanding, and the inertial modes of relativistic stars (of which the r-modes form 
a sub-class) have actually not been considered at all until very recently. In contrast, our understanding of rotating 
Newtonian stars has reached a relatively mature level and it is thus not surprising that most attempts to understand 
the r-mode instability and its potential astrophysical relevance have been in the context of Newtonian theory. 

Table [I] summarizes the differences between the low frequency modes of barotropic and nonbarotropic stars, and 
the the ways in which the relativistic inertial mode problem differs from the Newtonian problem (primarily because 
of the dragging of inertial frames) |3 ■ Following [j| , we use the term barotropic to describe a star for which the true 
equation of state describing both the background star and its perturbations is a prescribed one-parameter function 
p = p(e). In a nonbarotropic star the perturbations and background star obey different equations of state. The main 
cause of nonbarotropy in neutron stars is stratification via entropy or chemical composition gradients, with the latter 
being the most important for all but very hot (newly born) stars. 
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Nonbarotropic stars 


Barotropic stars 


Newtonian Theory 


infinite set [&] of r-modes for each [I, m] 
infinite set of g-modes 


a single r-mode for I = m 
infinite set of inertial modes 


General Relativity 


infinite set \8j of r-modes for each [I, m] 
infinite set of g-modes 
continuous spectrum? 


no pure r-modes 

infinite set of inertial modes 



TABLE I: A comparison of inertial mode results in Newtonian gravity and General Relativity. 

First, some notes on terminology. Perturbations of a spherical star can be decomposed into two classes depending on 
how the perturbed velocity transforms under parity (see 0). Following the standard relativistic terminology we will 
refer to perturbations that transform under parity like the scalar spherical harmonic Y" 1 as "polar" , while referring 
to those that transform opposite to Y™ as "axial" . This classification applies also to rotating stars, with the parity 
class of a mode being determined by its spherical limit along a sequence of rotating models |9j . 

Let us now characterise the rotationally restored modes; modes that have zero-frequency in the non-rotating (spher- 
ical) limit. We refer to modes that become purely axial in the spherical limit as r-modes, while modes that limit to 
a mixed polar/axial parity state are called inertial modes. Rotation breaks the degeneracy of these modes and gives 
them a finite frequency that is proportional to the star's angular velocity fl. We will also consider the g-modes, modes 
that are restored by the buoyancy associated with stratification. 

A rotating Newtonian barotrope possesses a vestigial set of I = m r-modes and a set of inertial modes, but no 
g-modes. A non-barotropic star, by contrast, possesses r-modes and g-modes. It is however worth noting that the 
low-frequency modes of a rapidly spinning neutron star may be similar to those of a barotropic model even though 
one would expect a realistic model to be stratified. If the Coriolis force dominates the buoyancy force, one would 
expect the low- frequency mode spectrum to be made up of inertial modes . Given that the g-modes of a "typical" 
neutron star model have frequencies below a few hundred Hz, it seems plausible that the low- frequency modes of 
millisecond pulsars will, in fact, be inertial modes. 

In Refs. 0| and we discussed the rotationally restored (inertial) modes of a slowly rotating relativistic star 
in some detail. One of the main results of this work was that these modes have a fundamentally different character 
in barotropic versus nonbarotropic stellar models. Firstly, the vestigial set of I — m r-modes that one finds in a 
Newtonian barotrope do not exist in a relativistic barotropic star pj. One is left with only inertial modes. This is 
particularly important because it is the I — m = 2 r-mode that is most likely to dominate the gravitational-wave 
driven instability. Secondly, it is possible to find r-modes of a relativistic nonbarotropic star at order £1 in a slow- 
rotation expansion. In the Newtonian case the r-modes are degenerate at order f2; one must undertake an order fl 2 
calculation to find the eigenfunctions. This degeneracy is partially split at order fl in the relativistic case, allowing 
one to compute the leading part of the mode eigenfunction. In the rest of this paper we will focus our attention on 
the relativistic problem for nonbarotropic stars, since the barotropic case is comparatively well understood and has 
been discussed elsewhere 

To lowest order in the slow rotation approximation, the r-modes of nonbarotropic stars are governed by an ordinary 
differential equation first derived by Kojima |12|. Consider the differential equation (for a complete derivation, see 
0) 
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(1) 



which determines the axial metric perturbations for a "pure" relativistic r-mode (h is directly related to Sgte)- In the 
equation v and A are coefficients of the unperturbed metric, and u~> is defined in terms of the relativistic frame-dragging 
cu as 



UJ 



Q - oj(r) 
ft 



(2) 



where Q is the (uniform) rotation rate of the star. Furthermore, we have assumed that (in the inertial frame) the 
mode depends on time as exp(icrt) and then introduced a convenient eigenvalue a as 
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Using the fact that 
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we have 



. . d h 9X , ^dh 

(a-u)< -—r - Anre^ip + e) — 
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— r + 8tt( P - 



e 2A h} + 167r(p + e)ae iA /i = (5) 



The above equation was first derived by Kojima an d our previous analysis Q shows that it can be used to 
determine r-modes of a nonbarotropic relativistic star. In order for the solution to satisfy the required regularity 
conditions at both the centre and at infinity, the eigenvalue a must be such that a — Co vanishes at some point in the 
spacetime 7]. As long as a — Q =/= inside the star, the problem is regular and one can readily solve it numerically. 
In our previous study we solved the problem for uniform density stars and found that the required eigenvalues were 
always such that the problem was non-singular |l3j |. There is of course no guarantee that the problem will remain 
regular for more realistic equations of state. Indeed, recent work by Kokkotas and Ruoff [Til fl5j and Yoshida 0] 
(see also 01) extends the analysis to more realistic equations of state, such as polytropes, and shows that the desired 
eigenvalue is then not generally such that a > uj s = Co(R). In other words, one is (at least for some stellar parameters) 
forced to consider a singular eigenvalue problem. Our analysis is intended to extend our results [jj for uniform density 
stars to more realistic equations of state. 

The existence of a singular eigenfunction problem has two consequences. Firstly, it has been argued that the 
singularity in Kojima's equation gives rise to a continuous spectrum of axial perturbations [l2L ITsj . One purpose 
of this paper is to explore the nature of the continuous spectrum and demonstrate that it is an artifact of the slow 
rotation approximation that may not be present in physical stars. It has also been suggested [l4l Il5lll^. Il7j that the 
r-modes may not exist in stars for which Kojima's equation is singular. Our aim in this paper is to argue that regular 
r-mode solutions will indeed exist in such stars. 

We will draw in part on studies of the oscillations of differentially rotating Newtonian stars, which also give rise to 
singular normal mode equations. Although the normal mode solutions are singular, the perturbation that one obtains 
by solving the initial value problem is non-singular [Hil I20I ] . The singularity is an artifact of assuming a normal mode 
time dependence. That consideration of the initial value problem regularizes the relativistic slow rotation inertial mode 
equations in the Cowling approximation has been demonstrated by Kojima and Hosonuma plj. We will show that the 
same is true for nonbarotropic stars when one includes the metric perturbations, despite the fact that the character 
of the singular normal mode solutions is different to that of the solutions found in the Cowling approximation. 

It is also plausible that the singular character of Kojima's equation represents a breakdown in the slow-rotation 
approximation. We will argue that one may regularize the singular normal mode equation in a physically well- 
motivated way by including higher order terms in the approximation. This issue, which develops further the work of 
Kojima and Hosonuma p^|. was addressed in some detail in an earlier preprint by two of us |23| . This paper is a slightly 
revised version of that original preprint. In particular, we incorporate recent improvements in our understanding of 
the nature of continuous spectra obtained from studies of differential rotation 0, [2(j . 

Other methods of regularization that have been discussed in the literature include the effect of gravitational radiation 
reaction, and coupling to higher order multipoles. For the nonbarotropic problem, gravitational radiation reaction 
alone is not sufficient to regularize the singular solutions found at lowest order in the slow rotation expansion 0, 
Il7^ . Coupling to higher order multipoles, however, does regularize the nonbarotropic problem both in the Cowling 
approximation 24] and when one includes the metric perturbations 22] . A brief note on the occurrence of continuous 
spectra and regularization in the barotropic problem is included in the Appendix. 



II. SINGULAR EIGENFUNCTIONS 



Let us consider Eqn. (JSJ in the case when a is such that we have 

a - Cb(r ) = (6) 

for ro < R, i.e. when the problem is singular at some point ro in the stellar fluid. Suppose we use a power series 
expansion to analyze the behaviour of the solutions to JSJ) in the vicinity of the singular point. Expanding in a; = r — ro, 
and assuming that all quantities that describe the unperturbed star are smooth, we use 

Co « tD(ro) + dj\x + Co-ix 2 ... where a — a>o(ro) = (7) 
P ~ po+pix + ... e w e + eix + ... (8) 
e A « e A ° [l + Aix-K..] M w M + M±x + ... (9) 

Next we introduce the Frobenius Ansatz 

00 

h = Y^ a n x n+(3 (10) 

n=0 
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in (J5J and find that we must have either j3 = or j3 = 1. This problem thus belongs to the class where the difference 
between the two values for (3 is an integer and we would not expect the two power series solutions to be independent. 
Indeed, further scrutiny of the problem reveals that we can only find one regular power series solution to our problem. 
This leads to an approximate solution 



/i rcg » a x [1 + a x x] 



where 



a x = 2tt{ Pq + e )e 2A ° ( r + ^ 



4a \ 



Cbi J 



(11) 



(12) 



In order to arrive at a second, linearly independent, solution we resort to the standard method of variation of 
parameters. Given one solution h\{r) to JSJ, a second solution can be obtained as 



h 2 = f(r)hi 

Introducing this combination in J5J it is straightforward to show that we must have 

f"hi= [47rre 2A + e)fti - 2h[] f 
(where a prime denotes a derivative with respect to r). In other words 



which integrates to 



f =^exp 



4irre 2X (p + e)dr 



(13) 
(14) 

(15) 
(16) 



Unfortunately, in our case we only know hi in the vicinity of the point ro . In order to proceed we therefore expand 
(|16fl in terms of x and then use hi = h r ° s . Then we need 



exp 



4nre 2X (p + e)dr 



1 + Eix 



where Ei — Anroe °(po + £o) and 



Putting the various pieces together we have 



/ no \ 2 _ 1 2qi 
\h rc &) ~ x 2 x 



dx 



C 



1 Ei - 2ai 
1 - 

2 



(17) 



(18) 



(19) 



with C an arbitrary normalisation constant. Integration then yields (recalling that x can take on both positive and 
negative values) 



f^-Ci- + (2a 1 -E 1 )ln\x\ 



(20) 



At the end of the day, we have arrived at a second solution to our problem (we discuss the consequence of taking 
In |a;| rather than In a; below). Near the point ro this solution can be written 



where 



h sin& ra b {l + b c x\n\x\+aix} 
167ra(p + £o)e 2Ao 



b c = 2ai — Ei 



CJl 



(21) 



(22) 
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(note that we need to keep the last term in 1)2 111 to work at an order that allows us to distinguish the leading order 
term of from the corresponding term in the singular solution). From this expression it is clear that, while the 
function /i slng is regular at r — r$ its derivative is singular at this point. 

In addition to this, one can show that it is not possible to find an overall solution to the problem (that satisfies 
the required boundary conditions at the centre and surface of the star) if one assumes that h oc /i reg in the vicinity 
of tq. Given this result we would seem to have two options: One option is to conclude that we must have a singular 
metric/velocity perturbation, and since this would be unphysical we must rule out the associated solution. If we 
take the im plic ations of this to the extreme, it could imply that no relativistic r-modes can exist for certain stellar 
parameters However, this conclusion is likely too extreme. It would be surprising if a small change in, say, 

the compactness of the star (the stiffness of the equation of state) could lead to such a drastic change in the star's 
physics (the disappearance of its r-modes). An alternative (and perhaps more reasonable) option is to assume that 
the appearance of a singular cigcnfunction signals a breakdown in our mathematical description of the problem rather 
than a radical change in the physics. Later in this paper we will show that the problem arises because of a breakdown 
in the slow rotation approximation. However, even in the slow rotation approximation, the perturbation is in fact 
completely regular; the presence of the singularity in Kojima's equation is simply a consequence of the assumption of 
normal mode time dependence. 

The normal mode equations for differentially rotating Newtonian stars exhibit mathematically identical singular 
behaviour for frequencies that lie within what we call the co- rotation band 0, |2(j • The eigenfunctions associated 
with this frequency band have singular derivatives that possess in general both a logarithmic singularity and a finite 
step in the first derivative at the singular point (see equations (45) - (49) of 0] and the accompanying discussion). 
The additional degree of freedom associated with the finite step in the derivative permits the existence of a continuous 
spectrum of solutions within this frequency band. At certain frequencies, the finite step in the first derivative vanishes; 
these frequencies are referred to as zero-step solutions and they possess a special character (see below, and the 
discussion at the end of Section 6.2 of |l9j')- 

The situation for Kojima's equation is identical: in general the singular eigenfunctions possess both a logarithmic 
singularity and a finite step in the first derivative, leading to a continuous spectrum of singular solutions. Just as in 
the differential rotation problem, there are certain frequencies for which the finite step in the first derivative vanishes. 
It can be shown that taking the logarithm of |a;| in the scries expansions in Eq. (20) and demanding continuity of 
the function at the singular point is equivalent to the matching procedure used in Section 6.2 of |19| to pick out the 
zero-step solutions from within the continuous spectrum of the differential rotation problem. Thus by using In \ x\ 
rather than In x in the analysis above we are picking out the zero-step solutions from the continuous spectrum. 

The perturbation is however determined by solution of the initial value problem rather than the normal mode 
problem. Analysis of the initial value problem for differentially rotating systems has shown that the perturbation 
associated with the continuous spectrum is not singular |2(j . By conducting a similar analysis of the time-dependent 
form of Kojima's equation, we have confirmed that the same is true for the relativistic r-modes. The singular solutions 
associated with the continuous spectrum are therefore physically relevant, and cannot be discounted. 

With this in mind, let us review the key characteristics of the differential rotation continuous spectrum and ask 
whether similar characteristics are manifested in the relativistic problem. Firstly, the continuous spectrum was found 
to possess a position-dependent frequency component; such behaviour has been observed in numerical time evolutions 
of the relativistic problem [T^j . Secondly, there were fixed frequency contributions from the endpoint frequencies of 
the continuous spectrum. Ruoff and Kokkotas |14| found such contributions in their simulations, but attributed them 
to the behaviour of the energy density at the surface of the star. We believe that they may instead be a hallmark of 
the continuous spectrum. The third characteristic of the continuous spectrum was a power law decay with time. In 
|14| there are two indications of this type of behaviour. The amplitudes of the endpoint frequencies were observed to 
die away as a power law. In addition, the authors noted that there appeared to be no contribution from the continuous 
spectrum at late times, suggesting again that it had died away. 

Consideration of the initial value problem for differential rotation also indicated a special role for the zero-step 
solutions |20j |. Again, the perturbations were found to be non-singular. For appropriate initial data the zero-step 
solutions were found to behave in much the same way as regular modes outside the co-rotation band, giving rise to 
a clear peak in the power spectrum at a fixed frequency and standing out from the rest of the continuous spectrum. 
The zero-step solutions behaved as modes within the continuous spectrum. Analysis of the time-dependent form of 
Kojima's equation indicates that the same will be true for the zero-step solutions to the relativistic problem. These 
solutions are therefore relevant. This contradicts statements in earlier works 0, fl5l liH Il7j that considered only 
the normal mode problem. The authors of these studies discounted these zero-step solutions within the continuous 
spectrum as being unphysical, and concluded that if r-modes entered the continuous spectrum they ceased to exist. In 
fact they do continue to exist as physically meaningful zero-step solutions, and should appear in time evolutions. For 
polytropic background models, Ruoff and Kokkotas |l4j| observe no contribution at the expected zero-step frequency 
when they initialise their simulations using arbitrary initial data. It would interesting to see whether these modes 
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could be excited using initial data more closely matched to the zero-step eigenfunction; the zero-step oscillations 
observed in |2("| were excited using initial data closely matched to the ei gen function rather than arbitrary initial data. 
For more realistic equations of state, however, the time evolutions of JlJ] do show clear peaks at fixed frequencies 
within the continuous spectrum. This suggests the presence of zero-step solutions. 

Before moving on we should make one comment on the nature of the continuous spectrum if one makes the 
Cowling approximation. In the Cowling approximation the continuous spectrum eigcnfunctions for the velocity 
perturbations are delta functions |2l| . Contrast this to the situation outlined above, where the velocity perturbations 
are proportional to the derivative of the metric perturbation, with a logarithmic singularity and (in general) a finite 
step at the singular point. The non-singular perturbation in the Cowling approximation (found by considering the 
initial value problem) exhibits a position dependent frequency component but no power law time dependence, no 
endpoint frequency contributions, and no zero-step solutions. The nature of the problem is changed dramatically by 
working in the Cowling approximation. 

We have argued above how the singular solutions of Kojima's equation give rise to non-singular perturbations when 
one considers the initial value problem. However, the main cause of confusion is a breakdown in the slow-rotation 
approximation. After all, Eq. (jSJl should really be written 



.d 2 h 



-l,Tre 2X (p + e)^ 
dr 



1(1 + 1) AM 
v 2 1 r + 8ir(p + e) 



e 2X h\ + 16Tr{p + e)ae^h = + O(n 2 ) (23) 



„2A 



From this we can immediately see that it is inconsistent to use the slow-rotation expansion when a — uj ~ 0(fl 2 ) or 
smaller. For the problem at hand this means that the assumptions used in the derivation of Eq. J^J) are not consistent 
in the vicinity of r$. Near this point we cannot discard the higher order terms while retaining the term proportional 
to a — uj since the latter becomes arbitrarily small. 

At first sight this may seem quite puzzling but similar situations are, in fact, common in problems involving 
fluid flows. In such problems, the singularity is usually regularized by introducing additional pieces of physics in a 
"boundary layer" near the point tq. A typical example of this, that has already been discussed in the context of the 
r-mode instability, is provided by the existence of a viscous boundary layer at the core-crust interface in a relatively 
cold neutron star (see Q] for an extensive discussion) . In that case the non- viscous Euler equations adequately describe 
the r-modc fluid motion well away from the crust boundary, while the viscous terms are crucial for an analysis of the 
region immediately below the crust. In our view, the relativistic r-mode problem leads to a similar situation: Well 
away from the point rg Eq. (JSJ leads to an accurate representation of the solution, but if we want to study the region 
near ro we need to include "higher order" terms in our analysis. 

Unfortunately, this means that it becomes very difficult to find a complete solution to the problem. The order fl 2 
perturbation equations for a relativistic star are rather complicated and have not yet been obtained completely. But 
for our present purposes, we can use partial results in this direction. Kojima and Hosonuma p2l | have shown that the 
next order in the slow- rotation expansion brings in a fourth order radial derivative of h in Eq. JSJl. Retaining only 
the principal part of the higher order problem we then find that © will be replaced by an equation of form 



uj 2 g(r)ar z ^ + (a 



d 4 h 



„ [d 2 h 



-\xre 2X (p + e)-^ 
dr 



1(1 



+ 1) 4M 

- — + 8n (P + e ) 



„2A; 



16Tr(p + e)ae 2X h = (24) 



where g(r) contains information about the stellar background — in particular the stratification of the star. Most 
importantly g(ro) 7^ and it is therefore clear that the problem is perfectly regular also near the point where a — Co = 0. 



III. A SUITABLY SIMPLE TOY PROBLEM 



Our main objective is to argue that one can in principle regularize the nonbarotropic r-mode problem. Ideally, 
we would like to find the mode-solutions without actually having to derive the relativistic perturbation equations to 
higher orders in the slow-rotation expansion. In other words, we are interested in a simple, practical approach to this 
kind of problem. 

As was shown in the previous section, the relativistic r-mode problem has (essentially) the following form 

fly"" + xy" + By = (25) 

in the vicinity of the point x = r — ro = (we use primes to indicate derivatives with respect to x). Both this toy 
problem, and the problem outlined below, retain the main character of Eq. (|24() but are sufficiently simple that we can 
solve them analytically. From standard perturbation theory, we know that this class of problems can be approached 
via matched asymptotic expansions. Typically, the outcome is that the singular equation (the equation obtained by 
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taking Q — > 0) leads to an accurate solution well away from x = 0, while the higher order term is required to regularize 
the solution near the origin. To illustrate this, and to motivate the method used to solve the r-mode problem in the 
next section, we consider the toy problem 

ey"" + x 2 y" - xy' + y = (26) 

where e is small in some suitable sense. 

Assuming a power series expansion in e we see that we first need to solve the singular equation, 

x 2 y" - xy' + y = 0. (27) 

The two solutions to this equation are y\ = x and y2 — xln\x\. In other words, the solutions to our toy problem 
are similar to the two (local) solutions we found for the relativistic r-mode problem in Sect. II. Hence, a method for 
solving our toy problem should be equally valid for the r-mode problem. 

Let us now suppose that we are interested in a global solution that satisfies boundary conditions y(l) — 1 and 
y'(l) = 0. Then we must have 

y{x) = x — xln \x\ (28) 

As was the case in Sect. II, this function is well behaved at the origin but its derivative diverges (cf. Fig. [IJ. Note 
that this solution also satisfies the boundary conditions y{— 1) = —1 and y'(— 1) = 0. 

Let us now consider the full fourth order equation i|26|) . It is straightforward to solve it using power series expansions. 
Inserting y = X^^o a nX n m Eq. I|26(l we find the recursion relation 

( n ~ 1 ) 2 Qn (0Q] 

" +4_ e(n + A)(n + 3)(n + 2)(n + l) ^ 

From this we see that we have four independent solutions. One of these, corresponding to a\ ^ truncates and leads 
to the solution y — a\x. As a result of the simple recursion relation, we can write the general solution to Eq. (|26[l in 
closed form: 

y(x) = a y Q (x) + a 1 y 1 (x) + a 2 y2{x) + a 3 y 3 (x) (30) 

with 



i=0 



Vi(x) = x (32) 

y v ' ^ (4z + 2) e' 1 v ; 

i—0 

where we have defined the symbol fc!!!! = k(k — 4)!!!! with k\\\\ = 1 for k < 0. 

We now want to find the specific solution to the fourth order problem which satisfies the boundary conditions 

y(-l) = -1 (35) 

2/(1)= 1 (36) 

?/(-!) = (37) 

y'( i) = (38) 

so that it agrees with the second order (singular) solution at the boundaries. It is straightforward to show that the 
required solution is ao = a 2 = and 

aa = , n f 1} m (39) 
2/3 CO -2/3(1) 

a 3 = tm r 1 m ( 4 °) 
y^ 1 ) -2/3(1) 
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xy' 3 (l)-y 3 (x) 

This solution is compared to the singular solution H28fl in Fig. ^ From the data shown in the figure one can conclude 
that the solution to Eq. 1|2()|) is well-described by the singular result (|28[) as long as we stay away from the immediate 
vicinity of x = 0. 




FIG. 1: We compare the singular solution of our simple toy problem to the complete solution (for e — 2 x 10~ 4 ) of the higher 
order (non-singular) equation. While the two functions agree well over the entire range (left panel), the derivatives obviously 
differ near the singular point at the origin (right panel) . The figures illustrate that the singular solution provides an acceptable 
approximation to the true solution well away from a region near the origin. 



IV. THE R-MODES OF NONBAROTROPIC RELATIVISTIC STARS 



The discussion in the previous two sections has crucial implications for our attempt to solve the relativistic r-mode 
problem for nonbarotropic stars. Clearly we can use our two solutions to Eq. JSJ to approximate the physical solution 
to the problem away from r = r$ even though one of these expansions is technically singular at rg. This provides us 
with the means to continue the numerical solution of Eq. JSJ across r — r$ , even though we will not be able to infer 
the exact form of the solution in a thin |25| "boundary layer" near this point. Should we require this information wc 
must carry the slow-rotation calculation to higher orders and solve a much more complicated problem. 

We thus propose the following strategy for solving the relativistic r-mode problem in nonbarotropic stars: First 
integrate the regular solution from the origin up to r = tq — <5, where 5 is suitably small. Then use the numerical 
solution to fix the two constants do and bo in the linear combination (cf. and (|21[l ) 

h = h ICS + h sins (42) 

Finally, this approximate solution is used to re-initiate numerical integration at ro + 5. This approach was first 
advocated by one of us in a set of circulated but unpublished notes |26) . and the idea was resurrected by Ruoff and 
Kokkotas [3. 

We have used the proposed strategy to calculate r- modes for a wide range of polytropic stellar models. Typical 
results are shown in Figure [21 (Shown also for comparison are the inertial mode frequencies of fully relativistic 
barotropes The inertial modes shown are those that limit to the I = m = 2 r-mode of the corresponding 

Newtonian barotropic model.) By comparing the obtained mode-eigenvalues a to the values for the relativistic 
framedragging at the centre and surface of the star (uj c and <D S , respectively), one can see that the r-mode problem 
is always regular for uniform density stars. As the equation of state becomes softer (n increases) the situation 
changes. For example, for n = 1 polytropes one must typically consider the singular problem in order to find the 
relativistic r-mode. This conclusion is in agreement with Kokkotas and Ruoff [T3as well as Yoshida [lg. It is worth 
emphasizing that the inertial mode problem for barotropic stars is never singular 0, ITlj unless one makes the Cowling 
approximation |2lLl27| . an approximation that is not in fact appropriate in the barotropic case. 
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FIG. 2: The r-mode eigenfrequencies a of relativistic nonbarotropic stars for n = (left panel) and n = 1 polytropes (right 
panel) . Also shown are the corresponding values of the relativistic framedragging at the centre ui c and surface d> 3 of the star. 
Whenever u) c < a < ci) s the problem is formally singular. As is clear from the data, the uniform density case (n = 0) is always 
regular while most of our n — 1 models are in the singular range. Also shown (as a dashed curve) are the eigenfrequencies for 
the axial-led inertial mode of a barotropic star that most resembles the leading Newtonian r-mode. Note that the inertial mode 
problem is never singular. 



Before discussing our results further we need to comment on a difference between our calculation and those in 
[T3. IT^ | . In these papers the authors consider polytropic equations of state of the form 



p = Ke 1 



-l/r. 



with p the pressure and e the energy density. Meanwhile, we are using 

and e = po + np 



T s 1 + 1 It, 



(43) 



(44) 



where po is the rest-mass density, in order to stay in line with the analysis of the inertial modes of barotropic stars 
[Tl| . This means that our numerical results cannot be directly compared to those in 0] . In order to verify that the 
results are consistent we have done some calculations using also (|43H . We then find that our results are in perfect 
agreement with those of Ruoff and Kokkotas. 

Our calculations thus support the numerical results of the previous studies. It is clear that, for more realistic 
equations of state one must consider the singular r-mode problem. Where we differ from both Ruoff and Kokkotas 
[l4| and Yoshida 0] is in the interpretation of the results in these cases. Yoshida only considers the regular problem, 
and tentatively argues that there may not exist any relativistic r-modes when the problem is singular. Similar 
conclusions are drawn by Ruoff and Kokkotas As we have already indicated, we disagree with these conclusions. 
Even in the slow rotation approximation the perturbations, obtained by solving the initial value problem, are non- 
singular. However, the root cause of the singular nature of the mode problem is a breakdown in the slow-rotation 
approximation. We believe that this problem would not arise if the calculation were taken to higher orders in f2 in the 
vicinity of the "singular" point (in analogy with boundary layer studies in problems involving viscous fluid flows). The 
physical problem is likely to be perfectly regular, but unless we extend the slow-rotation calculation to higher orders 
(or approach the problem in a way that avoids the slow-rotation expansion) we cannot solve the r-modc problem 
completely for nonbarotropic stars. However, we have shown how the r-mode eigenfrequencies can be estimated using 
only the solution to the singular mode problem, where they manifest themselves as zero-step solutions. 

The case in favor of our approach has been argued (we believe convincingly) in the previous sections. In addition, 
we can provide one further piece of evidence. In our previous study Q, it was pointed out that there is a striking 
similarity between the eigenfunctions of modes in barotropic and nonbarotropic stars. For example, the metric variable 
h(r) for an I = m = 2 r-mode of a nonbarotropic uniform density star was very similar to that of the axial-led hybrid 
mode corresponding to the Newtonian I — m — 2 r-mode. This is exactly what one would expect if the two represent 
a related physical mode-solution. We can now extend this comparison to the case of polytropic stars. The relevant 
data are shown in Figure We believe these data provide further support for the relevance of our nonbarotropic 
relativistic r-mode results. 
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FIG. 3: The eigenfunction for an r-mode of a relativistic nonbarotropic star is compared to the corresponding axial-led hybrid 
mode of a barotropic model. In each comparison, the equilibrium model is chosen to be the same: a relativistic polytrope of 
compactness M/R = 0.2 and polytropic index n = 0.5 (left panel) or n = 1.0 (right panel). The left panel shows a case in 
which the nonbarotropic mode is regular while the right panel shows a case in which the nonbarotropic mode is singular (the 
singular point is close to the surface at ro = 0.913i?.) The functions are all normalized so that h(R) — 1. 



CONCLUSIONS AND CAVEATS 



We have discussed the calculation of r-modes of relativistic nonbarotropic stars, shedding new light on a problem 
that has been associated with some confusion in the literature. We have shown how the seemingly singular problem 
can (in principle) be regularized, using standard ideas from boundary layer theory and viscous fluid flows, and how 
one can nonetheless estimate the eigenfrequencies of the desired r-modes from the singular mode problem. There are 
however issues that remain to be resolved, two of which merit particular comment. 

Kojima's equation admits a continuous spectrum of singular solutions whose collective perturbation is non-singular. 
The time-dependence of the collective perturbation is complicated, but includes a position dependent frequency 
contribution, and possible power law decay with time. At certain frequencies within the continuous spectrum one can 
find perturbations that behave like stable modes, whose manifestation is again non-singular (the zero-step solutions). 
We have argued in this paper that the underlying physical problem can be regularized by considering higher order 
rotational corrections. The effect of such regularization on the continuous spectrum and zero-step solutions is as yet 
unclear. If the zero-step solutions become regular normal modes then this would be physically interesting. The fate 
of the rest of the continuous spectrum is unknown; it may remain, vanish, or break up into discrete normal modes. 

The second issue is of relevance should we want to assess the astrophysical importance of the r-modes we have 
computed. In order to do this we need to estimate the timescale on which the mode grows due to gravitational 
wave emission |lll l'2<Sl . This calculation requires knowledge of the perturbed fluid velocity in order for the relevant 
canonical mode-energy to be evaluated. In the notation of Q, we need the variable U(r). We know from Eq. (4.21) 
in that 



[a-u))U 



-ah 



Clearly if we were to use our mode solution to Eq. the corresponding result for U would necessarily be singular, 
blowing up like l/(a — uj) at the singular point. In accordance with the arguments in Sections II and III above, it 
is easy to argue that the "physical" solution U(r) will be smoothed out by including higher order terms near the 
singular point and thus be regular at all points inside the star. However, solving this higher order problem is difficult. 
We can in principle avoid having to solve the higher order problem by solving instead the time-dependent initial value 
problem for the physical velocity perturbation. The physical velocity perturbation, just like the metric perturbation, 
will be non-singular. Unfortunately solution of the initial value problem is very difficult if one does not have a full 
analytic solution for the singular mode problem (see |20j where the same issues are discussed for the differential 
rotation problem). This may well mean that we cannot meaningfully estimate the gravitational radiation reaction 
timescale for the singular nonbarotropic modes discussed in this paper. 
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appropriate to describe the inertial modes of a barotropic star. The singular problem that arises in this case is thus 
an artifact associated with an unphysical assumption. 

In Ref. [27|, Ruoff, Stavridis and Kokkotas study barotropic inertial modes in the Cowling approximation. They 
expand the eigenfunctions in terms of spherical harmonics, which leads to a set of coupled equations that they truncate 
at some value l max of the angular parameter I. For a given l max , they find certain frequency bands for which the 
matrix problem cannot be inverted, and claim (correctly) that these frequency bands represent continuous spectra. 
When l m ax is increased, these continuous spectra are replaced by a discrete eigenfrequency solution at a frequency 
close to that previously occupied by the continuous spectrum. However, other continuous spectra now appear at 
different frequency bands. As l max increases these continuous spectrum bands grow in number and begin to span the 
full range of frequency space. 

The authors do not explain why this should be the case, but note that the continuous spectra may vanish when 
further higher order couplings are taken into account. In fact they should also vanish in the limit l m ax — * oo with 
only the lower order couplings that they consider. The continuous spectra that they observe for a given l max are 
the continuous spectra associated with the inertial modes of highest I. By including one more term in the coupling 
equations these solutions are regularized by coupling to higher order multipoles, hence the appearance of discrete 
mode frequencies. At the same time a new set of continuous spectra appear, associated with the unregularized 
inertial modes that have I = l ma x- As I increases there are more and more inertial modes [2^, hence the apparent 
proliferation of continuous spectra. In the limit l max — ► oo the inertial mode problem will be regular, and concerns 
that the frequency band may fill up with continuous spectra are unfounded even in the low coupling approximation. 
In this limit mode "disappearance" within continuous spectra is no longer an issue. In the limit of finite l max modes 
may however appear as zero-step solutions within the continuous spectra: this is a topic for further study. 



